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p ^. Abstract 

. We study the stability in finite times of tlie trajectories of interacting particles. Our aim 

[ is to show that in average and uniformly in the number of particles, two trajectories whose initial 

positions in phase space are close, remain close enough at later times. For potential less singular than 
the classical electrostatic kernel, we are able to prove such a result, for initial positions/velocities 
distributed according to the Gibbs equilibrium of the system. 



< 

: 1 Introduction 

The stability of solutions to a differential system of the type 

1^ 



^' 
O 
O 



is an obvious and important question. For times of order 1 and if F is regular enough (for instance 
uniformly Lipschitz), the answer is given quite simply by Gronwall lemma. For two solutions Z and 
to fll.ll) . one has 



\Z{t) - Z\t)\ < \Z{0) - Z\0)\ exp(t ||VF||l.o). (1.2) 

This inequality forms the basis of the classical Cauchy-Lipschitz theory for the well posedness of (11. ip . 
It does not depend on the dimension of the system (the norm chosen is then of course crucial). This is 
5^ I hence very convenient for the study of systems of interacting particles, which is our purpose here. 
Consider the system of equations 



vr = ENixn = ^ K{x[^ - xf) 



(1.3) 



where for simplicity all positions Xj^ belong to the torus T and all velocities belong to M . This 
system is obviously a particular case of (11.11) with Z = Z^ = {X^ , . . . , X^, V-^ , . . . , V^). 
The equivalent of (II. 2p reads in this case 

||Z^(t) - Z^'^(t)||i < ||Z^(0) - Z^'^(0)||i exp(t (1 + llVKlUoo)), (1.4) 
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where we define the norm on 11^^ x M.^^ 



1 ^ 



1=1 



This estimate is logically uniform in the number of particles A^. It is important in itself but also because 
it is a crucial tool to pass to the limit in the system of particles and derive the Vlasov-type equation 

dtf + v-VJ+(^F^,(^j^Jit,x,v)dv^^ ■V./ = 0, (1.5) 

for the 1-particle density f{t,x,v) in phase space, where * denotes the convolution. Hence estimates 
such as (11. 4p are at the heart of the derivation performed in [20], [lO], and [1] (we also refer to [23] . 
[23] and |26]). Note that the derivation of coUisional kinetic models (of Boltzmann type) involves quite 
different techniques, see [18] or [5]. 

Unfortunately, many cases of interest in physics deal with singular forces K ^ W^^^. Typical 
cases are K = — V0, a periodic force coming from a periodisation of the potential (/>R(a;) ~ C/|x|°~^, 
i.e. (f>{x) = (pmi^x) + g{x), where g{x) is an (at least) function on the torus T^. As the potential </> is 
defined up to a constant, we may also assume that its average is 0: J^g 0(x) dx = 0. The most important 
case is the electrostatic or gravitational interaction: a = 2, in dimension 3. 

Very little is known for these singular kernels, either from the point of view of the stability or of the 
derivation of Vlasov-type equations. Provided a < 1 and the initial configuration of particles are well 
distributed, the limit to Vlasov equation (11. 5p was proved in [Ti] . 

For systems without inertia, i.e. when the equations are simply 

{ = E^iXn = ^ K{X[^ - Xf) , (1.6) 

it seems to be easier to implement Gronwall-type inequalities. The derivation of the mean field limit is 
consequently known up to a < 2 (a < (i — 1 in dimension d), see [13] and also [T7] for a situation where 
the forces have a more complicated structure. In this setting the most important case is however found in 
dimension 2, for K = (corresponding to a = d — 1 = 1); the limit is the 2d incompressible Euler 

equation written in vorticity form. The derivation of the mean-field limit in this case was rigorously 
performed in [11] and [21], [22] . 

For differential equations like (II. ip in finite dimensions, it has long been known that well posedness 
and stability (for almost all initial data) can be achieved without using Gronwall-type estimates. The 
introduction of renormalized solutions by DiPerna-Lions in [9] gave well posedness for F G W^'^ with 
divF e L°°. 

This was extended to F G BV in the phase space situation in [3] and then in the general case in [1] 
(see also [15] for a slightly different approach). The exact case of the Poisson interaction was treated in 



This well posedness implies some stability as the flow has then some differentiability properties, see 
[2]. However the corresponding stability estimate is not quantitative and this kind of method does not 
seem to be able to provide uniform estimates in the number of particles (which gives the dimension of 
the system). We refer to [8] for a precise overall presentation of the well posedness and differentiability 
issues for Eq. (11. ip in finite dimension. 

More recently a new method to show well posedness for (II. ip has been introduced in [7]. Given a 
fixed shift 6, it consists in bounding quantities like 
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where Z is the flow associated to f ll.ip . i.e. Z(t, Z^) is a solution to fll.ll) satisfying Z{0, Z^) = 

A bound on such a quantity shows that for a.e. Z^ the two trajectories and Z° + S) 

remain at a distance of order |5| = |2'(0, Z*^) — ^(0, Z*^ + In this sense, this is an almost-everywhere 

version of the Gronwall inequality (11. 2p . 

It was shown in [7] that the quantity (11.71) remains bounded if F G W^'"^ for some p > 1. This was 

extended to W^'-^ and SBV in [TB] and even to H^l^ in the phase space setting in [B]. However in all 

those results the bound depends on the dimension of the space and is blowing-up as this dimension 

increases to oo. 

Therefore, our precise aim in this article is to prove a bound on a quantity like (II .yp for the system 
(II. 3p . uniformly in the number of particles. To our knowledge, this is the first quantitative stability 
estimate to be obtained for singular forces. 

Several new important issues occur when one tries to do that though. One of the most important is 
the reference measure which is chosen as this can imply different notions of almost everywhere as the 
dimension tends to +oo. In finite dimension, this refers to the Lebesgue measure and of course implies 
corresponding estimates for any measure which is absolutely continuous with respect to the Lebesgue 
measure. In infinite dimension, no such natural measure exists. This is due to the phenomenon of 
concentration of measures. Even in the case of finite but large dimensions, this is a problem to get 
quantitative estimates. Indeed even if two measures and V2 on 11^^ x are absolutely continuous 
with respect to each other or even more if 

dvx < C dh'2, 

then the constant C will in general depend on the dimension and go to +oo as increases. This 
means that a uniform quantitative estimate on the trajectories for some measure dvi{Z^) on the initial 
configuration will not in general imply a good estimate for another measure dv2{Z^). 

For each A^, the choice of the measure fi^ on U^^ x M.^^ is hence crucial to get a good estimate. One 
would naturally want to select a measure fi^ which is bounded, stable and invariant by the flow, just 
as the Lebesgue measure is stable and invariant by the flow of (ll.ip when F is divergence free. Let us 
denote by Z^{t) = {X^{t), (t)) or Z^{t, Zq) the vector of particles velocities and positions evolving 
through Eq. (ll.Sp till time t and depending on the initial configuration Zq . The system (II. 3p has an 
invariant which is the total energy 

Hn[Z^] = E^ + ^E'^d^'-^fl) (1-8) 

i i^j 
= Ek^n{V'') + EpotiX"") (1.9) 

To get an invariant measure /ijv, the simplest choice is to take a function of the total energy. Among 
those which are stable, the most natural is the Gibbs equilibrium 

rf^^(Z^) = ^e-^^-[^"lrfZ^ , (1.10) 

where dZ^ is Lebesgue measure on T'^^ x M.^^ , and 

S^(/3) = j e-^^-[^"l dZ^ , (1.11) 

is the normalization constant. Note that for a potential (p with a singularity at the origin, this makes 
sense only if limo^ = +oo, that is in the case of repulsive interactions. In the following, since we deal 
with measures which have a density with respect to the Lebesgue measure dZ^ , we will use the same 
notation for the measure fi^ and its density. 
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We study the quantity 



Qit) = I d^^{Z^) I .p^{Z^, 5) In ( 1 + ''^ f ' ) d5 , (1.12) 

where (5 at is a small parameter that will go slowly to zero when goes to infinity. 5 is a shift on the 
initial condition Zq 

Here iIjn ■ x M^^ ^ ^(T^^ x R^^) (where P(fi) denotes the set of probabilities on fi) is a 
probability valued function, so that it satisfies 



/ #^^(Zo^(5) = l, VZo^ G (T^ X ^^'^^ 



t/'at then then gives distribution of the shifts on the initial conditions, and the quantity Q(t) is averaged 
both on the initial conditions Zq and on the shifts 6. 

We now define the image measure of /i at by the shift distribution: 

fiN{Zo) = / I2n{Zo - 6)i/jn{Zo - 6,6) d6 . 



The crucial assumption will be that the image measure fiN remains "close" to the original Gibbs measure 
in the sense that 

3i^/3 > 0, such that VZq and A^, fiN{Zo) < K^fiNiZo) (1-13) 

with a constant independent of A^, but which may depend on /3. 
We will also use the weaker but very similar condition 

3K'f, > s.t. VAT, 3/3' < (3, s.t. VZq, ^n{Zo) < i^>^(Zo), (1.14) 

where /j]^ denotes the Gibbs measure with inverse temperature (3. The last condition is more general 
than the first, it allows to control the image measure by a Gibbs measure with bigger temperature. 

Remark 1.1 We mention that by definition of ftN, the measure 7ri\i{ZQ, Zq) = /ijv(^o)'^Af(^0) ~ ^o) 
is a transport from the measure fiN to fiN- In fact, ipN is (up to a translation of origin) what is usually 
called the desintegration of the measure tin with respect to its first projection fi^. However, we preferred 
our less standard presentation since we are more interested in hn and its shift ipN that in the precise 
image measure fiN- We mention the analogy to emphasize that our quantity Q is related to optimal 
transport. In fact 

Qn{0)= [ dfir^iZ^) [ dijN{Z^,6)\n(l+^^] >Wr,{fiN,fiN) 

J J5eT3JVxM3iv V On / 

where Wn is the transport for the cost In ^1 + ^ j . 

Conditions f ll.lSp and f ll.l4p are not explicit on Roughly speaking, they will be satisfied if ipN 
is chosen so that \Hn{Zo + 6) — Hn{Zo)\ < C if 6 E Supp \l/Ar(Zo, ■). This is reasonable since Hn is 
preserved by the dynamics, so that if the shift changes the energy too much, the original and shifted 
dynamics may be very different. But that simple and "reasonable" condition is not sufficient, we will 
really need a bound like (I1.13P on the image measure constructed with the shift. 

As the conditions are not explicit, we will provide in section [2] some examples of admissible shift 
distributions. The main result of the paper is a control on the growth on this quantity Q: 
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Theorem 1.2 Assume that (p > (pmin for some (pmm G ond that for some constant C , and a < 2 



CC ^jQ 'JQ 



OL+l ' 



Then taking 5^ = N'"^ for any e < 1 — a/3 and for all N > ^2^^y2 one has 

Ccp + Cacl 



Qit) <[l + {l + Kp) t + Q(0) , 

where Cj^ = e~t'^'"'", a is any exponent strictly larger than 2a/3, C constant (that can be made explicit), 
and Ca satisfies Ca < g^^^a ■ 

This theorem is not able to deal with the electrostatic interaction, a = 2; gravitational is of course 
out of question since repulsive potentials are needed. Note however that the electrostatic potential is 
just the critical case. The same result could be obtained in any dimension, with essentially the same 
proof. In dimension d, the condition would then he a < d — 1. The growth of Q is linear in time: note 
that this indeed corresponds to an average exponential in time divergence of the trajectories, analogous 
to (O). 

Roughly speaking, and provided that the average shift at time is of order 6n (or smaller), the 
theorem says that the average shift transported by the dynamics remains of order 6n- during the evolu- 
tion, and the control given is quite good. It is interesting to compare 6n to the minimal distance in the 
{X, V) space between two particles of a configuration, which is of order A^~3. Notice that it is always 
possible to choose 6n smaller than A^~3. Then if the order of magnitude of the initial shifts is smaller 
than A^~3, the theorem says it remains so at all time. This implies that the pairing of a particle in 
the configuration Z{t) with the closest one in the configuration Z^{t) is not very much affected by the 
dynamics: in this sense, there is not much "mixing". 

While the Gibbs equilibrium is the most natural choice for the measure fiN, others are possible. The 
proof would work for any measure fiN such that 

- fiN is invariant under the flow or ji^i^Z^ {t)) = H]^{Zq) for Z^ solution to ( 11. 3p 

- for all the /c-marginal of //at defined by = J ^^{Z^, Z^ ^)dZQ^ ^ satisfies: 

VZo, /i^(Zo)<C^ 

Obvious candidates are functions of the renormalized energy -^Hj^ but checking the bound on the 
marginals is not necessarily easy. 

Link with mean field limit. 

Finally, let us emphasize that this stability estimate does not answer the question of the mean field 
limit. Doing so would require to be able to deal with much more general measures /iat. More precisely 
if one can prove Th. 11.21 for a sequence of fiN and if in some reasonable sense 

- nilif 0, 

then the mean field limit is proved but only for the initial data Currently the Gibbs equilibrium 
corresponds to f^{x,v) = e~^^'"^ 

Unfortunately, we are not able do deals with more general measures /x^. The problem is that we 
need bounds on every k marginals and those are very difficult to obtain if we start from another measure 
than the Gibbs equilibrium. For instance, starting from fiN = g'^^ for some smooth profil g, we have 
the desired bounds at time 0, but do not know if fiN{t) satisfies them for any other time t > 0. 
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2 Some examples of admissible shift distributions. 
2.1 Shift on velocity variables 

In this section, we will be interested in shifts acting only the velocities. A first possibility is to take 
shifts independent of Zq and acting independently and identically on each velocity variable. Precisely, 
we are looking for shift distributions such as 

N 
i=\ 

where t/^ is a probability on symmetric with respect to the origin (■?/;(—?;) = if i\) has a density). 

We will not be able to deal with a general but will show that the hypothesis f ll.l4p is satisfied if 
■0 is Gaussian or has a compact support. This is stated precisely in the following Proposition: 



Proposition 2.1 Assume that is a Gaussian probability with variance a 



2 



then the hypothesis (I1.14p is satisfied with 



3 



f 1 \^ 

e 2^ 



/3'(iV) = /3 ( 1 - ^ ^ ) , and Kp = e-^<^-e 



2,2 „ „ 



4^- 



Assume that ip has a compact support with Supp-?/' C 5(0, 5m) (the ball of radius 6m, center 0). Then 
f ll.l4p is satisfied for N > with 

I3\N) = (3[l- --f , andKp = e ^ 62'^'^ 



Remark that for a < |, the average velocity fluctuation given by such shift is larger than the smallest 
5n we can choose. 

Proof of the proposition 12.11 

In the Gaussian case, we have 

The 3A^ integrals may be performed independently, using the ID calculation: 



We finally get 

ZN 

3N 
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Epot(Zo) 



with (3\N) = - j^j^). Using 



1 + N/(3a 
we get 

We define B^^x the normahzation integral restricted to the position variables: 

Bn,xW) = J e-^^''°*(^i'-'^^) dXi... dXN 
From the Jensen inequality applied to the function x i— > x^'^^, we obtain for l3' < l3 

Then, using the bounds of Lemma [3. ![ we have Bn,x{P'Y' ^ 1; this implies 
This proves the desired inequality with 

In the case of ip with compact support, we follow the same sketch. To do this, we will need a bound 

on 



Using the symmetry of iIj, and the inequality cosh(x) < e^2 (valid for a: G M), we may bound that term 
by 

cosh \^--j=j e" 2iv dip{5v) < J e ■2n''-2n dip{5y) < e 
and as in the previous case, we get 

1 -b( HN(Zo)-^Ev.i„(Zo)] 



B 



N 



From this point, following exactly the same step as in the case of gaussian ip, we prove that the hypothesis 
(11.141) is satisfied with the announced constant □ 

By making the shift depend on the velocity Vq = (Vi, . . . Vat), one may essentially remove the con- 
dition on the size of the norm of the shift. More precisely we limit ourselves to shifts 6 = (0, . . . , 0, S') 
with 5' G R^^, giving Zo + 6= (Xq, Vo + 6') with Xq = (Xi, ...,Xn). Now define 

iPn{Zo, 6) = {Ul,6s^=o) ^{\6'\) G{\Vo\) S,Vo-s'+\sr=o, (2-1) 
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where = ^-^ is the usual euchdian norm and 5... is the corresponding Dirac mass on the hypersurface 
of equation 2 Vq ■ 5' + = | Vq + — |Vop = which is precisely the sphere of radius |Vo|. 
We need the function ipN to satisfy 

1= [ijN{Zo,6)d6 = G{\Vo\) I 

J ^2 Vb-<5' + |<5'|2=0 

this is always possible with the right choice of G as the integral 

J2Vo-5' + \5'\^=0 

depends only on |Vo| by the rotational symmetry of the sphere. 

Condition (11.131) is automatically true since, as //at depends only on Xq and |Vb| and |Vo — 5'| = |Vb|, 

f^NiZo) = j iin{Zq - 5) iPn{Zo -5,S) = Hn{Zo) j i^NiZ^, -5) = fiN{Zo), 

because iPn{Zo — S,6) = ipN^Zo, —6). One could wish to impose additionally that \6'\i < 6n, so that 
Q{0) is of order 1. Since 

\S'\, = l^y\6l\<^\6l 

i ^ 

it is enough to impose that ^ has support in [0, A^^/^ 5^]. 
As a conclusion, we proved 

Proposition 2.2 For any measure onM.^^ , there exists a function G{\Vq\) , s.t. the probability density 
ipN defined by ( 12. ip satisfies ( 11.13^ . 

One could try to generalize this example by letting Vo + 6' and Vq to be on close but different energy 
spheres. For example by posing 

MZo,S) = {nl,Ss.=o) ^m)G,{\Vo\) Ii2Vo-s+s^\<v 

Provided that t] in not too large, the computations are essentially the same and one has essentially to 
make sure that 

G^,{\Vo\±v)<GG,{\Vo\). 
2.2 Shifts in position variables 

One could try to implement the same idea for shifts in position variables. Many problems arise however 
since the potential energy is not a nicely regular function of the positions. 

If one tries to consider shift distributions ipNi^) that do not depend on Zq, then the limitation on 
|5| is quite drastic. In fact this example essentially works if only a fixed (independent of N) number of 
coordinates 6i are not 0. 

Trying to generalize the second example by imposing that Zq + S and Zq are on the same energy 
surface also faces several problems. The main one is that the equation of the energy surface is not 
anymore symmetric between the shift S and the shift —6. 

The only solution would be to write the equation only on the tangent plane, i.e. something like 

= '^m)SvHiZoy5=o\VH{Zo)\. 
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This is now nicely symmetric but poses other difficulties. For instance one would need to make sure 
that H{Zo + 6) < H{Zq) + C . Expanding H , one would formally find a condition of the type 



Unfortunately V^-ff(Zo) is singular and in particular unbounded. It is only bounded in average which 
would force us to remove the initial conditions around which the energy has a singular behavior. 

Although this procedure could in principle be carried out successfully, we do not wish to enter here 
into such technical computations. This essentially limits us to present explicit shift distributions acting 
only on velocity variables. 

However, let us point out that the evolution of the particles system strongly mixes positions and 
velocities. Obviously, if we start from two initial conditions Zq and Z^ with same positions and different 
velocities, we get at time t > configurations with different positions. 

So if we really need a shift distribution that acts also on the positions at the origin t = 0, a strategy 
may be to start at t = — r > 0, using a shift distribution acting only on velocities at this time, and then 
to let evolves the particles till time 0. The original shift distribution transported by the flows is now 
a shift distribution acting on position and speed. Using the Theorem (11. 2p . we get that the average of 
these evolved shifts (in a weak sense since we we are taking a logarithmic mean) is at most of order 5ni 
provided the average of the original shifts was also of order at most 5n- Therefore by removing a set of 
vanishing measure of initial conditions, one obtains a shift distribution in positions and velocities that 
satisfles all the requirements. 

3 Some useful bounds for the 6A^ dimensional ji^ 

We shall make use of the following lemmas. Before stating them, we introduce a notation for the 
projection of //at on the position space. 



V'H{Z,)\\5\'<C. 




(3.1) 



We also denote the k- marginal of vjsi by 



Lemma 3.1 For all N, we have: 




(3.2) 



We will also need the following estimate on the kth marginal of /iat deflned by 




Bn 



1 



I 



e 



Uz 



■N-k 



(3.3) 



and more precisely on the k marginal z/^ of the projection un on the position variables. 




(3.4) 
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Proof of lemma [3. II To compute the integral defining Bn fll.lip . we may separate the integration in 
from the integration in . In the V'^ variable, we have to integrate a product of 3A^ independent 
real gaussians of variance (3~^. We obtain (27r//3)^^/^. 

In the variable, we use Jensen inequality by the convexity of exponential to get: 

r, \ ~3N/2 p 



which gives the first inequality (We used that has zero average). To obtain the second bound, it 
suffices to use the inequality Epot{ZQ) > ^(pmin- □ 

Proof of lemma 13.21 The proof follows the one introduced in [19] for the Lame-Enden equation. 
As the measure /iat factorizes in position and speed, we may write 

Bn,x 

Neglecting the terms in interaction energy involving (at least) one of the first k particles, we obtain 

Bn,x J 



B 



N,X 



the term ^^^j^ being there because Epot{Xk) = (1/^) Yli^j 'Pil^i ~ ^jl) ^ positions. So we need an 
estimate on terms of the kind 



e{k) = j e-^^'^^^^^'^'Ux' 



We can relate this term to configurations with k + 1 particles. First use Jensen inequality as the 
exponential is convex to get 

Q{k)=j e-/5^^-*(^'=)dX^ 

exp||-y" (^fi^Ep,,{X'^) + ^J2^{\X,-Xk+i\?j dxk+?j dX^ 

<Q{k + l). 

Since, 0(iV) = Bn,x, we iterate this inequality and get Q{N — k) < Bn,x- Putting this in (13. 5p . we get 
which is the result needed. □ 
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4 Proof of Theorem 11.21 

During the course of the demonstration, C will denote a constant (independent of and which 
value may change from line to line. 

From now on, we shall omit the superscript in the notation , as there will be no ambiguity. 
We have to estimate the derivative of Q{t). Differentiating directly, one obtains 



^^^w-y \ Sn+\\Z-Z^\\^ 



+ 



6n+\\Z-Z% 



where Z^ = {X^, V^) = Z{t, Zq + 6). 

Note that the first term is obviously bounded by 1 and hence 

dnu,^. , f, fv, ,,wE.\E,mX,-X,)-KiX^-X^))\ 



5n + \\Z-Z^\\i 
We define for a integer L that will be fixed later 

Ci{Zo,t) = {j 7^ i s.t. \Xi{t) — Xj{t)\ is among the L smallest \Xi — X^l} (4.1) 

That is for each i, Ci collects the indices of particles which are closest to particle i at time t, following 
the flow. It also depends on the initial condition Zq. We define similarly Cf. 
Accordingly, we decompose dQ/dt as follows 

j^Qit) <C + S^ + Si + S2, 



with 



S, = j dfi^iZo) j d6MZo,6) j^Yl E \K{X,-X^)\, 

S( = Jdf^^iZo) Jd6iJNiZo,6)j^J2 E \KiX^-X^)l ^4 2) 



r r 1 \K{Xi- Xj) - K{x° - 

4.1 Bound on Si 

Recalling the bounds on i^' = — V0 from Th. 11.21 one simply begins with a discrete Holder inequality 
for any 7 < 3 

1— a/7 

I r r 1 _ _ \ 

^1 < 



dfiN{Zo) I d6ijN{Zo,6) ^ 1 j 
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We first use the fact that the integral of V'jv in 5 is equal to 1 to get 



We then perform the change of variable from Zq to Z, using the inverse fiow of Z{t, Zq) which preserves 
the measure fi^; one finds 



As the second marginal of fi^ is bounded by ci by Lemma [3 .2^ this implies 



X,|7 



EEn^ 



X,|7 - (3-7)5 



N 



Hence 



^1< 



(3 - 7)5iv 



To conclude, simply note that by definition |Cf | = |Cj| = L and so for any a > 



where Ca satisfies Ca < 



c 



3a— 2a ' 



4.2 Bound on Sl 

Using Fubini and the change of variable Zq y-^ Zq + 5, and the image measure vn-, we may rewrite 



And from the hypothesis fll.l3p . we may bound that that exactly as the previous one. The only difference 
is that a constant K\ appears and that we shall use /3'{N) instead of /3. We get 

4.3 Bound on 5*9 



5n \N 



By using the assumption on the second derivative of in Th. II. 2[ one first bounds 

|K(X, - X,) - K{Xf - Xf)\ < C - Xf I + \X, - X^l) + -xs^a^i 

Therefore defining the following matrix, with Ia the characteristic function of the set A: 

12 



one has ^ ^ 



Consequently, if we use the classical matrix inequality ||Mx||i < sup^ ( ^ ■ |Mjj|)) ||x||i, 



5*2 < y" duN^Zo) j d5iljN{ZQ,5) max^Mj^. 



As before the terms in M containing \Xf — are the equivalent of the ones with \Xi — Xj\, thanks to 
(11.131) . Hence one has to bound 



S2 < C (5*2 + S2J, 

with 



5*2= I dfiN{Zo) I d6tpiy{Zo,6) max X] _ ^ 
S2 = J djjNiZo) J d5ip]\}{Zo,5) max 



a+l ' 



j s.t. i^Cj 



Since nothing depends on 5 now, one may integrate ipNiZojS) in 6 with value 1. Moreover changing 
variable from Zo to Z (we recall the flow is measure preserving), one simply finds 



Si = I dfiN^Z) max ^ \x - \ 
Si = / dfir,{Z) max Yl 



j S.t. i^Cj 



Let us now carefully bound each of these terms. 

The term 

We use 



/r+oo 
f{X) dfiM = P{f{X) >l)dl, (4.4) 



where P is the probability with respect to the measure fi^ on U^^ x M'^^. We have to evaluate now 
expressions like 



P max 



-iy, — ^--—]>i]=pi3zs.t.-\y- — ^_|>/| . (4.5) 



\jic, I * ^' / / \ Kite, I ' -^1 



To bound this probability, we will need the following lemma: 

Lemma 4.1 Given xi, . . . , > and I > 0; given {uk)k=i such that ^^^^ Uk = I- IfYl^=i ^« > ^^^^ 
3k n], . . ,ik / Xi^ > luk , Vr = 1, . . . , /c. 

Proof: Let the Xj be sorted Xi > X2 > . . . > x„, and suppose the conclusion is not true. Then we have 

Xi < lUi, . . . , X„ < lUn ■ 

Thus Xli < ^ Ei Uk = l- □ 
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We apply the lemma to Eq. 14.5^ with Uk = i'^)'^ < z/ < 1 to be determined. Cn is chosen 
such that J2k=i'^k = 1- Using Riemann sums, we see that lim+ooCAr = 1, and that Cn > 1,VA^ > 1. 
Hence we get 




k = l ^ I « Jrl / 

N~L / . ^ . l/(a+l) / ^ \ A/3 



fc=l 
7V-L 



fc=i 



where for simplicity we have introduced the parameter A = 3(1 — v)/{a + 1). 

To estimate the probability Pi^ki once the particle i is chosen, we have a constraint on the position 
of k particles, which have to be close enough to particle i, plus constraints on the position of L distinct 
particles, from the definition of Cj. This event concern k + L + 1 particles, and to estimate it, we will 
use an estimate of its volume P^**^ in the configuration space T3(fc+L+i) ^ ]]^3(fc+L+i)_ thus involves the 

{k + L + 1) marginal of which is bounded by c^^^^^ by Lemma [3^21 
This leads to the following estimates 

Moreover, using a simplified version of Binet formula (See [25] ) 

n! = v^n^'^^e-^+il^r , for some 9 G (0, 1) 
the binomial coefficient may be bound by: 

N\ , , 1 On Op «iV-P f N , 

fjP = = l27rl 2pl2iV 12p 12(JV-p) ( (4 fil 

^ p\{N-p)\ ^ ' pP{N - p)^~P \p{N - p) J ^ ' 



1 1 



N y fNY p A"-" fNeV 



< (2,)-,e™^_J + <^_J . (4,7) 

And we do not forget that since = C^~^ we may use the same inequality with p replaced by — p. 
Inserting this in the above inequality, we get: 



3(fc+L) 



k + L \vl \N 



< CcpN 



a + 1 



N J 



N \ ul 
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where = c^"'*'^'*^^, and A = e*^""*"^)/^ is a numerical constant. Now taking u close enough to (precisely 
z/ < ^^), one has A > 1 and then we take as well L > {X — 1)^^ (recall that L has yet to be fixed) ; 
hence 

3(k+L) 

.'Ac'„\ "+i 

Pi,k < Ccfs k 

If we sum on /c, we get: 



ul 



fc=0 ^ ^ fc=0 ^ ^ 

■A{L+1) 3(Z.+1) 

provided I > Iq = -. If we take moreover L > p > p{a + l)/3 for some p > 0, we get a simpler 

bound: 



k=0 

Remark that the conditions on L depend only on a and A (which depends only on a). In particular, 
those conditions are independent of the parameter /3. Thus, we have 



< ^4^"°^' (z//)~P for / > I'q (4.11) 



Integrating this quantity in /, one obtains for M such that Cc^ ^ {vM)~^ = 1 that 

/■oo /"Af /"oo 

/ P{--->l)dl = / P{--->l)dl+ P{--->l)dl 

Jo Jo Jai 

< M + Cc,3 



< M + Ccp 



Ac'pY r dl 
Jm 



V ) p —\ 



p- 1 - z/ ^ 

for p > 23^- And finally, we get that 

for some constant C, provided that L is large enough. 

Since every calculation have been performed, we see that a possible for v \^ v = in which case 
the condition on L is exactly 

6 

L > 

- 2-a 

With that choice of u, we get 

^2^ < ^ , 
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The 5*1 term: Through the same type of computations, we are led to evaluate expressions like 



< f p(3,;.,....W.^C.l^-i-^>.^|^(±)") (4.13) 

< E^(3j;^i,---,^^/j-^ai^.-Xd<^^(^-j j (4.14) 



fc=l 

N 



(4.15) 



k=l 



Since the sum is performed on the particles i such that j ^ Ci, we cannot choose L + k particles close 
to j as for 5*2. But, we have nevertheless to choose k particles close enough to j, a probability that will 
give a good bound if A; > L. If k < L, once a particle i close to j such that j ^ Ci is chosen, one knows 
that there exist L other particles close to j. This will be enough to bound the probability. 

In the second case {k < L), we pick up a particle j {N possibilities), at least another particle i (since 
k > 1) and then we have to choose L other particles closer to i than j is. Since \Xj — Xi\ has to be less 
than (A;/iV)(i-'^)/("+i), 

L L . . 3L/(q+1) / l \ AL 

fc=0 k=0 ^ ^ ^ ^ 

^ /^^X (A-l)L /^^/ X 3L/(a+l) 



fc=0 



2 ^E^A'M^ h? (^-1^) 



< CcJ^^'l^) f^) . (4.18) 



If L < V A^, we may use 



L 



.N 

6 



(A-l)L 



as soon as L > In that case, we may use this bound in last inequality and obtain: 

j:pl < cci (^) 

fc=0 ^ ^ 



/ ^ 3L/(a+l) 



if L < TiV. 

In the first case k > L, we pick up the particle j, and then choose k particles v close to j. We obtain 
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as previously 

N~l _N_ ^ ^ ^ 3fc/(a+l) / \k 



— r_i_i i,_r_i_i \ / 



k=L+l k=L+l 

N 



N 



^ N {A-l)fc ^^^/ X 3fc/{«+l) 



k=L+l ^ ^ ^ ^ 



fc=L+l 

^ / \ 3fc/(a+l) 



^^^^ 3L/(a+l) 



fc=L+l 

where we again restricted ourselves to (A — 1)L > 1 and assume / > /q • Putting the two sum together, 
we get the bound 

N-l /^^/^ 3L/(a+l) 



k=i ^ 



ul 

It remains to integrate in I. Doing exactly as for the term, and choosing the same t^, we will get 

The only difference is that it will require p > and thus L > 

4.4 Conclusion of the proof 

Putting all together, we finally we may bound 

dQ . ^ ( lV'^^ , c 

I — I J 



with C/3 = e"^*^""" and a > where L is subject to the restrictions (with the choice v = which 

means that A = 4t-^ ) 
1+2,0. ' 

L>^, , L<Vn. 
2 — a 

It is possible only if > that case it is clear that L should be chosen as small as 

possible and from the constraint that means 

L 



2- a 



With this choice, one has 



dt - ^2-a'^^ 6j,N'-^ ^ 2-a 
Now if one takes 5n = N~^, we can get a uniform bound in A^, only if 
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If this is true, we get 

— ^ < IH -. 

at 2 — a 

with Ca < 3^^2a which is the result given by Theorem 11.21 
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